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Abstract 

Roads and vehicular traffic are among the most pervasive of threats to biodiversity because they fragmenting habitat, 
increasing mortality and opening up new areas for the exploitation of natural resources. However, the number of vehicles 
on roads is increasing rapidly and this is likely to continue into the future, putting increased pressure on wildlife 
populations. Consequently, a major challenge is the planning of road networks to accommodate increased numbers of 
vehicles, while minimising impacts on wildlife. Nonetheless, we currently have few principles for guiding decisions on road 
network planning to reduce impacts on wildlife in real landscapes. We addressed this issue by developing an approach for 
quantifying the impact on wildlife mortality of two alternative mechanisms for accommodating growth in vehicle numbers: 
(1) increasing the number of roads, and (2) increasing traffic volumes on existing roads. We applied this approach to a koala 
(Phascolarctos cinereus) population in eastern Australia and quantified the relative impact of each strategy on mortality. We 
show that, in most cases, accommodating growth in traffic through increases in volumes on existing roads has a lower 
impact than building new roads. An exception is where the existing road network has very low road density, but very high 
traffic volumes on each road. These findings have important implications for how we design road networks to reduce their 
impacts on biodiversity. 
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Introduction 

With landscapes becoming increasingly dominated by humans, 
the impact of roads on wildlife populations is growing rapidly 
[1,2]. The direct effects of roads are wide ranging and include the 
destruction and modification of habitat [2], the modification of 
animal behaviour [3], the fragmentation of habitat by the 
formation of barriers [4,5] and vehicle collisions [6] . Roads also 
indirectly affect wildlife populations by increasing human access to 
previously inaccessible areas and changing land use patterns [2,7]. 
Consequendy, understanding the impact and the conservation 
management implications of roads is a high priority for a wide 
range of species of conservation concern and for concerned 
scientists [8]. 

A particularly important impact associated with roads is 
elevated mortality rates from vehicle collisions. This has been 
shown to have substantial impacts on wildlife populations [6,9,10]. 
Fahrig et al. [1 1] show that the density of anurans (frogs) decreases 
and the proportion of dead anurans increases with traffic intensity 
on roads near Ottawa, Canada. They conclude that high traffic 
volumes increase mortality enough to significandy reduce popu- 



lation densities. Similarly, Jones [12] demonstrates the extinction 
of a population of eastern quolls (Dasyurus vivemnus) in Tasmania, 
Australia following a road upgrade. She links this to an increase in 
vehicle collision mortality due to higher vehicle speeds. Therefore, 
management of the mortality effects of roads is a critical 
consideration for the conservation of the ever-increasing propor- 
tion of biodiversity that occurs in close proximity to human 
setdements. 

A major challenge for mitigating the impact of road mortality 
on wildlife is that, in almost all parts of the world, the number of 
vehicles on roads is increasing rapidly, with consequent increases 
in ecological impacts [13,14]. Despite the importance of managing 
vehicle growth in a way that minimises impacts on wildlife 
populations, we currendy have little understanding about the 
consequences of alternative road network design strategies for 
limiting the impact of vehicles on wildlife populations. This is 
because the focus to date has been predominantly on either only 
quantifying impacts [15] and/or evaluating mitigation measures 
such as road overpasses/ underpasses and fencing on existing road 
networks only [16-18]. Evaluations of the implications of 
alternative road network designs are much rarer [19—21]. It is 
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therefore important that we address this gap by placing more 
emphasis on evaluating alternative road network designs if we are 
to make informed decisions about future road network design 
strategies. 

In terms of planning road networks, there are essentially two key 
ways in which the growth in vehicle numbers can be accommo- 
dated: (1) by upgrading existing roads to carry higher volumes of 
traffic and/ or (2) by increasing the number (or density) of roads in 
the network. However, these two strategies result in different 
spatial distributions of traffic and therefore contribute to wildlife 
mortality rates via two different processes. If we increase the 
density of roads, but keep the traffic volume on each road 
constant, higher mortality rates arise because the probability that 
an animal moving around the landscape will cross a road increases 
[6]. On the other hand, if we increase traffic volume on each road, 
but keep the density of roads constant, higher mortality rates arise 
because the probability that an animal crossing a road is hit by a 
vehicle increases [6]. Although, under either strategy, mortality 
rates increase, it is far from intuitive which one results in the lowest 
increase in mortality, because that will depend upon the relative 
impact on mortality of two quite different processes. For a 
decision-maker faced with a choice between the two strategies, it is 
therefore critically important to understanding which strategy is 
best, and under which circumstances. 

Friar et al. [19] present one of the few examples where 
alternative road network designs are evaluated with respect to road 
density and the placement of roads in relation to habitat. 
However, they do not consider how their results vary with traffic 
volume. On the other hand, van Langevelde and Jaarsma [20] do 
explore road network design strategies that modify the spatial 
distribution of traffic volumes, but only on existing roads. Jaeger et 
al. [21] consider the impact of alternative road network designs, 
including one where traffic volume is concentrated along a single 
road, but they only explore this in simple artificial landscapes. 
Therefore, an explicit evaluation of road network design strategies 
that modify road density versus strategies that modify traffic 
volumes in real landscapes is a research priority. 

Here we address this issue using a spatially-explicit simulation 
model to quantify the relative impact of changes in road density 
and/ or traffic volume on mortality rates for a koala (Phascolarctos 
cinereus) population in eastern Australia. For koalas, mortality on 
roads can form a large component of overall mortality rates in 
many areas and it is considered to be one of the key threatening 
processes for this species [22-26]. We characterise model outputs 
using a statistical approximation and ask whether general 
principles emerge about the relative benefits of accommodating 
more traffic by increasing the density of roads (i.e., more roads) 
versus increasing traffic volume on existing roads (i.e., larger 
roads). We show that, under most circumstances, it is preferable to 
accommodate a greater number of vehicles by increasing the 
capacity of existing roads, rather than building new roads. 

Materials and Methods 

We used an existing model of koala movement for the study 
area [27] and no ethics permits or permissions were required to 
undertake the study. 

Study Species 

The koala is a folivorous and arboreal marsupial restricted to 
the eucalypt forests and woodlands of eastern and south-eastern 
Australia. Across its range, the koala feeds on a wide variety of tree 
species, predominandy from the genera Eucalyptus and Corymbia, 
but in any particular area, they show preferences for just a few 



species [28-32]. Koala habitat generally consists of the preferred 
food tree species in any area, although other factors, such as tree 
size and water availability, can also contribute to habitat quality 
[32-35]. Koalas occupy reasonably well defined home ranges and, 
although largely solitary, both male and female home ranges can 
overlap [30,36,37]. The key threats to the species are habitat loss 
and fragmentation, urbanisation, dog attacks, vehicle collisions, 
disease, bushfire and climate change [22,23,25,26,38-41]. Move- 
ments between trees, particularly in fragmented peri-urban areas, 
are usually made along the ground and this includes movements 
across roads. 

Study Area 

Our study area was located within the Port Stephens Local 
Government Area, New South Wales, Australia, approximately 
150 km north of Sydney (Figure 1). This area contains one of the 
most significant koala populations in New South Wales [42]. 
However, since European settlement, habitat loss and fragmenta- 
tion has occurred in the area due to agriculture, urbanisation and 
sand mining, bringing important threats to the koala population in 
the region [40,41,43]. Among these threats, vehicle collision 
mortality on roads is considered to be one of the most significant in 
Port Stephens [44] . 

Simulation Model 

To model the impact of roads on koala mortality, we combined 
a simulation model of koala movement with a model of the risk of 
mortality when crossing a road. Below we describe the model in 
detail and then describe how we applied it to our study area to 
quantify the relative impact of increases in road density versus 
increases in traffic volume. 

Movement Model. We used a spatially-explicit habitat 
selection and movement model to simulate koala movements in 
the study area [27]. Movement was simulated on a raster grid, 
with a 50 mx50 m cell size, representing the distribution of koala 
habitat in Port Stephens. The spatial distribution of habitat was 
derived from an existing koala habitat model and vegetation maps 
for the area [27,31,33]. The habitat model was developed based 
on information on koala tree species preferences estimated from 
field surveys of koala faecal pellets and this was then combined 
with detailed vegetation maps of the study area to arrive at the 
final habitat map (see Appendix C in Rhodes et al. [27] for a full 
description of the habitat mapping procedures). Each raster cell 
was classified as either: (1) primary/secondary habitat; (2) marginal 
habitat; (3) other vegetation not classified as koala habitat; (4) 
cleared; or (5) water bodies (Figure 1). Rhodes et al. [27] separate 
primary and secondary habitat, but we combined these two 
habitat categories here to reduce the number of habitat classes 
[41]. For a detailed description of the habitat categories see 
Lunney et al. [33]. The spatial habitat data is available from the 
Dryad Digital Repository: http://clx.doi.org/10.5061/dryad. 
3n4h2. 

Based on the movement model described in Rhodes et al. [27], 
we assumed that the probability of moving from location a to 
location b, Pr(a to b), in a given time period, in a landscape 
consisting of k= 1, m discrete grid cells, each of a defined 
habitat type j= 1, n, is 

<p(a,b) ± Wj(b)I(bJ) 

Pr(atob) = — -, (1) 

A ]T q>(a,c k ) £ Wj(c k )I(c k J) 

k=l ,; = 1 
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Figure 1. The Port Stephens Local Government Area. Map shows the study area's location in Australia, the estimated distribution of koala 

habitat, and the estimated average daily traffic volume (axle-pairs day -1 ) on major roads. 

doi:10.1371/journal.pone.0091093.g001 



where <j>{a,b) is a function defining the probability of moving from 
location a to location b independent of habitat (similarly for 
$(a,Ck)); Wj(b) is the relative preference for habitat type j at 
location b (similarly for Wj{Ck)); A is the area of each grid cell; Ct is 
the location of the centre of grid cell k; and I(bJ) is an indicator 
function which equals 1 if the habitat at location b is of type j and 
equals 0 otherwise (similarly for I{ctj))- The numerator defines the 
preference for moving from location a to b and the denominator is 
a normalisation constant that standardises this preference to a 
probability. The movement model represented by equation (1) is 
essentially a biased random-walk model with movement biased 
toward preferred habitat. It achieves this by modelling an 
underlying random-walk described by the function which 
defines the probability of movement to any location in the 
landscape in the absence of habitat selection, that is then modified 
by habitat preference through the function w,(). Note here that 
habitat preference is assumed to be dependent on the location of 
the habitat (see below). 

The habitat independent movement probability function, 
<j>(a,b), was defined as 



distributed. The habitat preference parameters, Wj(b), were 
defined as 



Wj(b)=exp(a.j + f}r bk ) 



(3) 



<p(a,b) = 



Aexp(-Xr a , 
2nr„h 



(2) 



where X is the scale parameter for the negative-exponential 
distribution and r ab is the distance between location a and location 
b. This assumes that, in the absence of habitat selection, the 
probability distribution of movement steps is negative-exponential 



where Oj is the preference for habitat j; r hh is the distance from 
location b to the centre of the animal's home range; and the 
parameter (3 defines how habitat preference varies with distance 
from the home range centre [27]. Negative values for ft imply a 
tendency to move back towards the centre of the home range and 
therefore introduce a form of home range behaviour. This model 
therefore enables the simulation of movement that accounts for 
both habitat selection and movement behaviour within a home 
range. For a full description of the movement model, habitat 
models and parameterisation see Rhodes et al. [27]. 

Road and Traffic Volume Data. Data on average daily 
traffic volumes (axle-pairs day ) from traffic recording stations in 
Port Stephens between 1995 and 2001 (New South Wales Roads 
and Traffic Authority unpublished data, Port Stephens Council 
unpublished data) were used to estimate average daily traffic 
volumes on all the major roads in the study area. These traffic 
volumes were then mapped spatially (Figure 1). By linking traffic 
volumes to the spatial location of roads, spatial variation in traffic 
volumes were explicitly defined for incorporation into the model. 
The traffic volume spatial data is available from the Dryad Digital 
Repository: http://dx.doi.org/10.5061/dryad. 3n4h2. 

For some traffic recording stations we also had data on hourly 
road traffic volumes from 2001 (New South Wales Roads and 
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Maritime Services unpublished data) and this revealed a dramatic 
difference in traffic volumes between day and night. Koalas are 
more active during the night-time than during the day [36,45], so 
the vast majority of road crossings will occur during the night. 
Therefore, we were interested in quantifying the proportion of 
traffic volumes that occur between 1800 h and 0600 h. We 
characterised this by fitting a beta distribution, by maximum- 
likelihood, to the proportion of daily traffic volumes occurring in 
each hour between 1800 h and 0600 h for the roads where we had 
hourly volume data. This provided estimates for the rate 
parameter, r, and scale parameter, s, for the beta distribution, 
beta(r, s) [46]. We estimated r= 1.32 and 5 = 87.16 and 
subsequently used this distribution, combined with the average 
daily traffic volumes, to draw random values for hourly traffic 
volumes on each road during the times when koalas are most likely 
to cross roads. 

Mortality Risk Model. To estimate the risk of mortality 
when a koala crosses a road, we used a simple model of the 
probability of being hit by a car. We assumed that the number of 
cars passing along a road per unit of time is Poisson distributed; a 
reasonable assumption for night-time vehicle volumes [47]. When 
a koala crosses the lane of a road, a gap between vehicles greater 
than the amount of time it takes the koala to cross the part of the 
road traversed by vehicles is required for a successful crossing. 
Therefore, assuming that koalas arrive randomly at a road and 
cross immediately, the probability of surviving a crossing is 

Pr{survival} = Pr{gap>Ar} = exp(— yAf), (4) 

where At is the time taken to cross the part of the road traversed by 
vehicles and y is the rate parameter for the Poisson distribution, 
representing the traffic volume. For a road with n lanes, with total 
two-way traffic volume, y, and assuming traffic volume is divided 
equally between the lanes, i.e. y/n on each lane, then the 
probability of surviving a crossing of the entire road [48] is 

Pr{survival}= [expf-Af^l =exp(— yAt). (5) 



This general derivation holds for any n-laned road, provided 
vehicle volumes are split evenly among lanes [47]. 

We further assumed that all vehicles on each lane travel along 
the same part of the road and that koalas cross perpendicular to 
the flow of traffic. In this case, the time taken to cross the path 
traversed by vehicles is 



where W is the vehicle width; / is the koala head to tail length; and 
v is the velocity at which koalas cross [47,48]. Hels and Buchwald 
[9] take a similar approach, but include the possibility that road 
crossings occur at different angles and that the killing width of the 
car is only a proportion of the vehicle width. The subjects of their 
study were amphibians, for which individuals are often only hit if 
they are under the wheels of a vehicle. We made the reasonable 
assumption, given the size of a koala, that passing anywhere under 
a vehicle would result in a fatal collision. For simplicity, and in the 
absence of data on the distribution of crossing angles, we also 
assumed that all crossings are made perpendicular to the flow of 
traffic. 

Movement Model Parameter Estimates and 
Uncertainty. Model parameters for the movement model for 



males and females were estimated from a koala radio-tracking data 
set for the Tomago Sandbeds region of Port Stephens [27] . This 
provided estimates of: habitat preference for marginal habitat, 
oc marg , other vegetation /mining revegetation, a,,^, cleared land, 
%ciem-, the negative-exponential scale parameter, X, and the 
parameter determining the influence of the distance to the home 
range centre, for each sex (Table 1). Habitat preference 
parameters were all estimated relative to primary/secondary 
habitat and water bodies were assumed unavailable. Head to tail 
lengths, /, were estimated as 0.66 m for females and 0.70 m for 
males using data for 69 adult female and 48 adult male koalas 
from South East Queensland (Table 1 , Queensland Department of 
Environment and Heritage Protection unpublished data). The 
velocity, v, at which koalas cross roads is uncertain and, to our 
knowledge, no empirical data currently exists for this. In the 
absence of such information, we assumed that this parameter 
would be somewhere in the range 5000-15000 m h" 1 (Table 1). 
These values are based on personal observations (D. Lunney) that 
koala movement velocities across roads would be at least as fast as 
a typical human walking speed of around 5000 m h -1 , but could 
be higher. We assumed that the average width of vehicles, W, was 
2 m (Table 1). 

We also characterised uncertainty in the estimates of the 
parameters a marg , cr.„ tkm d c y Mn X, [}/,„ I and v. We did this for the 
parameters a marg> oc alher , c/. clear , X and /} kr by assuming that they were 
distributed multivariate normal with expected values and vari- 
ance-covariance matrix estimated from the log-likelihood of the 
movement model [27], Koala head to tail lengths were also 
assumed to be normally distributed with means and variances 
derived from the sample means and standard errors of the head to 
tail length data [49]. Although we used sampling distributions, 
rather than Bayesian posterior distributions, to describe parameter 
uncertainty, when using uninformative priors, sampling distribu- 
tions will tend to approximate the Bayesian posterior distributions 
[50,51]. Therefore, these distributions were deemed to be 
adequate approximations for parameter uncertainties given the 
data underlying their estimation. The high degree of uncertainty 
in the velocity at which koalas cross roads, v, was characterised by 
considering three separate values for this parameter: 5000, 10000 
and 15000 m h" 1 . 

Simulations. For each simulation run, 500 male and 500 
female koalas were located randomly on the landscape and their 
movements were simulated for 365 days in 24-hour time steps. We 
assumed that their locations at the start of the simulations were the 
centres of their home ranges (i.e., they started from their home 
range centres). For each road crossed during a koala movement, 
we determined whether a mortality event occurred or not. This 
was achieved by first drawing a random variable from the beta 
distribution describing the probability distribution of hourly 
proportions of traffic volumes between 1800 h and 0600 h. This 
random proportion was then multiplied by the average daily 
volume for the road crossed to obtain an hourly traffic volume at 
the time of the crossing, thus assuming that the time of crossing 
was chosen randomly between 1800 h and 0600 h. Assuming that 
one axle-pair is one vehicle, the probability of surviving the road 
crossing was calculated from equation (5) using the randomly 
drawn traffic volume. Whether the individual survived the crossing 
was based on a draw from a Bernoulli distribution with probability 
equal to the probability of surviving the road crossing. An 
individual that did not survive a road crossing was recorded as 
such, but movement simulations were continued for the entire 365 
days so that an estimate of the individual's home range could 
subsequently be calculated. 
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Table 1. Parameter values used in the simulations. 





Parameter 


Symbol 


Baseline Value 


Standard Error/Range 


Females 


Habitat preference for marginal habitat 


rj -marg 


-0.210 


0.057 


Habitat preference for other vegetation 


a other 


-0.474 


0.086 


Habitat preference for cleared land 


C-clear 


-0.717 


0.118 


Negative-exponential scale parameter 


A 


4.77x10~ 3 rrT 1 


0.1360x10~ 3 


Distance to the home range centre parameter 


far 


-3.77x10~ 3 rrT 1 


0.1466x10~ 3 


Head to tail length 


1 


0.659 m 


0.008 


Movement velocity 


V 


10000 m rT 1 


5000-15000 


Males 


Habitat preference for marginal habitat 


V-marg 


-0.262 


0.073 


Habitat preference for other vegetation 


rj -other 


-0.396 


0.120 


Habitat preference for cleared land 


V-clear 


-0.373 


0.175 


Negative-exponential scale parameter 


X 


2.52x10~ 3 


0.0928 x10~ 3 


Distance to the home range centre parameter 


ftr 


-2.52x10~ 3 


0.1014x10~ 3 


Head to tail length 


1 


0.700 m 


0.012 


Movement velocity 


V 


10000 m rT 1 


5000-15000 



The habitat preference parameter for primary/secondary habitat was fixed at zero, so all habitat preference parameters are relative to primary/secondary habitat. cc marg 
= habitat preference parameter for marginal habitat, a otrier = habitat preference parameter for other vegetation/mining revegetation and a. ciear = habitat preference 
parameter for cleared land. 
doi:1 0.1 371 /journal.pone.0091 093.t001 



The simulated daily locations for each individual were then used 
to construct a 95% fixed kernel home range using a smoothing 
parameter, h, equal to the resolution of the landscape of 50x50 m 
[52] . Within the estimated home range, the road density (i.e., the 
proportion of grid cells in the home range that contained a road), 
X dms , and the mean traffic volume (i.e., then mean road traffic 
volume in grid cells containing roads), X mh were recorded. 

Statistical Analysis 

We used the simulation outputs to develop logistic regression 
models of the risk of mortality due to vehicle collision as a function 
of road density and traffic volume [53]. The response variable for 
these models was the binary mortality/ survival data for each 
individual, with X dens , X mi and a X dms by X m i interaction as 
explanatory variables. We excluded those individuals whose home 
ranges did not contain any roads. If there were no roads in an 
individual's home range, the mortality risk was always zero, so we 
were interested only in developing models for estimating the 
mortality risk, conditional on at least one road being present in an 
individual's home range. The regression models took the form 



In 



I -p. 



z <P + Vdens^dens + <P m [X vo l + <Pdens x vol^dens^vol > (?) 



where p is the probability of mortality and <j>, <Pj ens> (p vo i and 
fdenswol are the regression coefficients. Collinearity between mean 
traffic volumes and road densities was low based on Pearson's 
correlation coefficients (p<0.1) and therefore collinearity was not 
considered an issue for the regression models. 

To validate the models, we first ran simulations for each of the 
three alternative values of movement velocity, v, with baseline 
parameter values otherwise, for each sex (Table 1). Logistic 
regression models were fitted to each of the six resulting simulated 
datasets and we tested the model fits using Hosmer-Lemeshow 
deciles of risk and Pearson # 2 global goodness-of-fit tests [53]. The 



/(-values for the Pearson / tests were calculated from a normal 
approximation of the statistic's distribution [53-55]. 

We then used a bootstrap approach to estimate the expected 
values and standard deviations of the logistic regression param- 
eters so as to capture the uncertainty that arises from uncertainty 
in the simulation model parameters [56,57]. For each sex, and for 
each of the three values of v (5000, 10000 and 15000 m h" 1 ), we 
chose 100 random values oia marg , a. ot/ler , a rfean X, pi, r and / from the 
distributions describing the uncertainty in these parameter values. 
For each combination of parameter values, we ran a set of 
simulations and fitted the logistic regression model to the data, as 
for the baseline parameter case. This was repeated for each value 
of v and, for each of these, the mean and standard deviation of the 
estimated regression model parameters from the 100 replicates 
were calculated. This bootstrap approach provided, for each value 
of v, an estimate of the expected value and standard deviation of 
the regression coefficients [58] . The standard deviations primarily 
reflect parameter uncertainty, which is propagated through the 
random draws from the parameter distributions, rather than 
simulation error. Although some simulation error will be present 
in these estimates, the large number of individuals simulated for 
each parameter combination (500 of each sex) means that 
simulation error is likely to be relatively small. Efron and 
Tibshirani [58] recommend that between 50 and 200 bootstrap 
replicates are usually required to reliably estimate standard 
deviations. Therefore, the 100 replicates we used in this study 
were sufficient to obtain reasonable estimates of the parameter 
means and standard deviations. The bootstrap expected values of 
the regression coefficients were then used to make predictions 
about the probability of mortality due to vehicle collision mortality 
for values of X^ between 0.002 and 0.2 and X m! between 30 and 
30000 axle-pairs day . 

Finally, we aimed to quantify the impact of changes in X dens and 
X„i on vehicle collision mortality risk. To do this, we calculated the 
sensitivity and the elasticity of ln(p/l— p) to changes in road 



PLOS ONE | www.plosone.org 



5 



March 2014 | Volume 9 | Issue 3 | e91093 



Designing Road Networks for Wildlife Conservation 



density and traffic volume. Sensitivity with respect to road density, 
sj m „ and traffic volume, s m (, were calculated, from equation (7), as 



Sdens — 



df (Xj ens ,X ro l) 



and 



oyol : 



df(Xci em ,X vo i) 
dX m i 



z 9 dens + 9 dens x vol^wl, (8) 



= <Pvol + 9 dens x volXdens, (9) 



where f(X dem ,X m l) =<j> + 9dens X dens + 9vol X ™l + 9 dens x vol X densX m l 

is the linear predictor in equation (7). 

However, the values of road density and traffic volume were on 
quite different scales, making the relative interpretation of 
sensitivities difficult. As an alternative approach, elasticities 
provide a means of comparing the effect of proportional changes 
in variables, thus making them comparable. More specifically, 
elasticities provide a measure of the proportional change in one 
variable in response to a proportional change in another [59]. 
However, because the proportion of a logit has no intrinsic 
meaning, we modified the formulae for calculating elasticities so 
that we obtained the absolute change in ln(/)/ 1 —p) resulting from 
a proportional change in road density or traffic volume. In so 
doing, the elasticities with respect to road density, ej e „„ and traffic 
volume, e mh were calculated as 



t'di'l 



and 



: Xde, 



df (Xdens^Xvol) 
dX d ens 



' Xdens\9dens 9 dens x vol vol ). (io) 



e m i = X K 



df {Xdens,X vo l) 



(IX, 



vol 



-- X vo i((p m} + (p dens x vo iXdei 



(11) 



To compare relative elasticities, we used the ratio e dms /e m i and 
calculated this for X dms between 0.002 and 0.2 and X ml between 30 
and 30000 axle-pairs day , using the bootstrap expected values of 
the regression coefficients. To investigate robustness to parameter 
uncertainty, we also calculated the sensitivity, elasticity and the 
ratio of elasticities for each bootstrap replicate at typical mean 
values of X icns =0.04 and X ml =9500 axle-pairs day -1 . These 



values were then summarised by their expected values and 
standard deviations. 

Results 

Statistical Model Adequacy 

There was no evidence of a significant lack of fit for any of the 
logistic regression models fitted to the simulated data for the 
baseline parameters values based on either the Hosmer-Lemeshow 
deciles of risk or the Pearson y'f global goodness-of-fit tests 
(p>0.05). Therefore, the logistic regression models were consid- 
ered adequate descriptions of the relationship between mortality 
risk, road density and traffic volume. 

Regression Coefficients 

The expected values of the logistic regression coefficients 
showed that the probability of a mortality event on a road was 
positively related to road density and traffic volume, as expected 
(Table 2). There was considerable variation in estimates due to the 
propagation of uncertainty in simulation model parameters, but a 
high proportion of parameter draws resulted in positive slopes. 
However, the effect of road density was positive more often than 
the effect of traffic volume; coefficients for roads density were 
positive between 93% and 100% of the time, while coefficients for 
traffic volume were positive between 72% and 83% of the time. 
The coefficients for the interaction term was also positive, 
indicating that the impact of changes in road density was greater 
when traffic volumes were high than when they were low and/ or 
the impact of changes on traffic volume was greater when road 
density was high than when it was low (Table 2). Uncertainty in 
the simulation model parameters fed through to relatively high 
standard deviations (Table 2), but again, the proportion of 
parameter combinations that resulted in a positive interaction 
term was high (between 91% and 100%). Overall, positive 
coefficients occurred more often for males than for females. 

Predictions 

Predictions based on the regression models showed that males 
had higher annual mortality rates than females but, for both sexes, 
even low road densities were capable of causing high probabilities 
of mortality unless traffic volumes were very low (Figure 2). The 
major roads in Port Stephens have relatively high traffic volumes 
(the lowest recorded traffic volume on the roads used in the study 
was 1600 axle-pairs day 1 and the highest was over 30000 axle 



Table 2. Expected values and standard deviations (in parentheses) of the logistic regression coefficients for each sex and each 
movement velocity, v. 





V 


Intercept 


Density 


Volume 


Interaction 






(Xdens) 


(Xvol) 


(Xdens x X vo |) 


Females 


5000 


-2.56 (0.71) 


33.69 (17.85) 


4.25 x10~ 5 (6.89 x10~ 5 ) 


3.71 x10~ 3 (2.35x10~ 3 ) 


10000 


-2.97 (0.75) 


29.76 (17.98) 


5.12x10~ 5 (6.13x10~ 5 ) 


2.27x10~ 3 (1.91 x10~ 3 ) 


15000 


-3.09 (0.78) 


23.70 (15.53) 


3.28 x10~ 5 (6.65 x!0~ 5 ) 


2.28 x10~ 3 (1.54x10~ 3 ) 


Males 


5000 


-2.20 (0.55) 


43.18 (20.55) 


3.99x10~ 5 (5.14x10~ 5 ) 


6.67 x10~ 3 (3.17x10~ 3 ) 


10000 


-2.62 (0.52) 


35.98 (15.77) 


3.69 x10~ 5 (4.42 x10~ 5 ) 


4.59 x10~ 3 (1.90x10~ 3 ) 


15000 


-2.90 (0.58) 


33.05 (18.42) 


4.23 x10~ 5 (4.41 x10~ 5 ) 


3.53x10~ 3 (1.75x10~ 3 ) 



Values are the sample means and standard deviations of the parameter estimates for the 100 bootstrap replicates. 
doi:1 0.1 371 /journal.pone.0091 093.t002 
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Figure 2. Contour plot of the predicted annual probability of 
mortality. Values are shown as a function of mean traffic volume, X m i 
(axle-pairs day -1 ) and road density, X dens >0 (proportion of grid cells 
containing a road) for: (A) females and (B) males. Annual probabilities of 
mortality were calculated from the bootstrap expected values of the 
regression coefficients with v= 10000 m h _1 (Table 2). 
doi:10.1371/journal.pone.0091093.g002 

pairs day ) indicating that, even areas with low road densities 
may impose high rates of mortality on koalas. Mortality rates 
increased with road density and traffic volume, but generally more 
rapidly with road density than traffic volume (Figure 2). 

At typical mean values for road density (0.04) and traffic volume 
(9500 axle-pairs day - ), the sensitivities and elasticities were higher 
for males than females and declined as movement velocity, v, 
increased (Table 3). However the ratio of elasticities was similar 
across sexes and v and consistently showed a greater elasticity to 
proportional changes in road density than traffic volume. The 



Table 3. Expected values and standard deviations (in 
parentheses) of the sensitivities and elasticities of the logit 
probability of mortality with respect to road density, s dens and 
£dens and traffic volume, s vol and e vot and the ratio, e dens /e vo i, 
for each sex and each movement velocity, v. 







V 


^dens 


s vol 


®dens 




^dens^vol 


Females 


5000 


68.93 (12.22) 


1.91 x10 4 
(0.56 x10~ 4 ) 


2.76 (0.49) 


1.81 (0.52) 


1 .60 (0.37) 


10000 


51.37 (10.44) 


1.42x10~ 4 
(0.39 x10~ 4 ) 


2.05 (0.42) 


1.35 (0.37) 


1.61 (0.46) 


15000 


45.36 (10.35) 


1.24x10~ 4 
(0.34 x10~ 4 ) 


1.81 (0.41) 


1.18 (0.32) 


1.69 (0.78) 


Males 


5000 


106.57 (18.99) 3.07x10 4 
(0.95 x10~ 4 ) 


4.26 (0.76) 


2.91 (0.91) 


1.55 (0.38) 


10000 


79.57 (13.04) 


2.20x10 4 
(0.50 x10~ 4 ) 


3.18 (0.52) 


2.09 (0.48) 


1 .57 (0.34) 


15000 


66.62 (11.35) 


1.84x10~ 4 
(0.43 x10~ 4 ) 


2.66 (0.45) 


1.74 (0.41) 


1.60 (0.43) 



Sensitivities and elasticities were calculated at a typical mean road density of 
0.04 and a typical mean traffic volume of 9500 axle-pairs day -1 . Values are the 
sample means and standard deviations of the sensitivities and elasticities for 
the 100 bootstrap replicates. 
doi:1 0.1 371 /journal.pone.0091 093.t003 



elasticity with respect to road density was approximately 50% 
higher than for traffic volume. Bootstrap estimates of the 
sensitivities and elasticities were greater than zero 100% of the 
time, indicating strong support for an increase in mortality due to 
both increasing road density and traffic volume. There was slightly 
less strong support for the elasticity to road density being greater 
than the elasticity to traffic volume, but the ratio of elasticities was 
still greater than one for between 94% and 100% of parameter 
combinations. 

The elasticity ratio across the range of different values of road 
density and traffic volumes showed that, for most of the road 
density and traffic volume state-space, mortality was more elastic 
to a proportional change in road density, Xj m „ than the same 
proportional change in traffic volume, X ro/ (Figure 3). The 
exception was when road density was very low and traffic volume 
was very high, in which case there was a greater sensitivity to 
traffic volume. This general pattern held across all three values of 
movement velocity, v, and for both sexes, but especially for males. 

Discussion 

As the world becomes increasingly urbanised and human 
population sizes increase, identifying strategies to accommodate 
these changes while limiting impacts on biodiversity is critical 
[60,61]. Increases in road density and traffic volumes are typical 
features associated with greater urbanisation and human popula- 
tion growth, and are therefore fundamental considerations in 
conserving biodiversity [62] . A key planning decision that needs to 
be made for reducing the impact of roads on wildlife is whether we 
should accommodate increased traffic loads by increasing traffic 
on existing roads, by increasing the density of roads, or through a 
combination of both [63]. We have shown that, by using a simple 
model of animal movement and road mortality, key insights can be 
gained about the relative impacts of increases in road density 
versus increases in traffic volume on existing roads. In the vast 
majority of cases, we found that increasing road density elevated 
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Figure 3. Plot of the regions where the annual probability of 
mortality is more elastic to road density and where the annual 
probability of mortality is more elastic to traffic volume. These 
regions are shown as a function of mean traffic volume, X vol (axle-pairs 
day -1 ) and road density, X dens >0 (proportion of grid cells containing a 
road) for: (A) females and (B) males. Elasticities were calculated from the 
bootstrap expected values of the regression coefficients with 
v= 10000 m fT 1 (Table 2). 
doi:1 0.1 371 /journal.pone.0091 093.g003 

mortality rates more rapidly than did increasing traffic volume on 
existing roads. Decisions about where and how to build additional 
road capacity are governed by a range of factors, such as the 
spatial distribution of capacity requirements and implications for 
congestion [63]. However, our study indicates that strategies that 
focus on the creation of new roads are likely to be more harmful to 
wildlife than those that build capacity within an existing network. 

Previous studies have demonstrated that both road density and 
traffic volume can have substantial impacts on mortality rates and 



population dynamics [6], but we still know comparatively little 
about the relative impact of changes in road density versus traffic 
volume. Our most important new insight here is that the effect of 
road density on mortality is commonly higher than the effect of 
traffic volume. The opposite was only true when road density was 
very low and traffic volume very high. Therefore, it is only in cases 
where the existing road network is characterised by very few high 
capacity roads that increases in network capacity by building new 
roads is likely to have the least impact on koala mortality. The 
model also indicated that male koalas were more susceptible to 
road traffic mortality than females and the range of conditions 
under which building new roads was the better strategy was even 
more limited for males than females. The reason for this is that 
males generally have larger home ranges and move greater 
distances than females, particularly during the breeding season, 
with the result that males tend to cross roads more frequently than 
females. This is consistent with empirical evidence on differences 
between sexes in vehicle collision mortality rates in koalas [22,23]. 
The higher movement rates of male koalas also makes them more 
susceptible to the effect of increased road densities versus increased 
traffic volumes on existing roads. Therefore, in the case of males, 
this further reduces the range of conditions under which it is 
preferable to increase road network capacity by building new 
roads. The more general implications of this are that, for mobile 
species, upgrading existing roads is even more likely to be the 
better strategy than it is for less mobile species. 

Despite uncertainty in the parameter estimates for the 
regression models, the sensitivities and elasticities of mortality 
with respect to road density and traffic volume were greater than 
zero for all simulation model parameter combinations. Although 
we are slightly more uncertain about the relative effects of 
increases in road density versus increases in traffic volume on 
existing roads, mortality is still more elastic to road density than 
traffic volume for the vast majority of the simulation model 
parameter combinations. Therefore, even after accounting for 
parameter uncertainty, the conclusion that accommodating 
increased traffic through higher volumes on existing roads has a 
lower impact on mortality than accommodating increased traffic 
through the building of new roads is relatively robust. 

Although we focus on a single species here, the design of road 
networks will most commonly need to consider impacts on 
multiple species. Therefore, the extent to which our results can be 
generalised to other species is an important consideration. For 
example, species differ in their avoidance responses to roads, their 
movement speeds and their visibility to motorists and these factors 
may influence the relative impact of road density versus traffic 
volume [47,64,65]. We have shown that our conclusions are 
robust to movement speed, but we did not consider how avoidance 
behaviour or visibility to drivers could influence our results. 

Behavioural responses to roads have been observed for a range 
of species, particularly avoidance behaviour [2]. Avoidance of 
roads varies among species and can depend on a range of factors 
such as, the size of roads, traffic volumes and road noise [66-69] . 
Further, avoidance can be an important factor determining road 
mortality and may therefore modify the relative effects of road 
density and traffic volume on mortality [64] . Jaeger et al. [2 1] 
relate road avoidance behaviour to traffic density in an artificial 
landscape and show that the impact of roads on persistence is 
lower when traffic volumes are concentrated along a single road 
than when they are not. Although they only consider an artificial 
landscape, their results are consistent with ours across a range of 
traffic volumes and avoidance probabilities. This indicates that our 
results may be robust to assumptions about avoidance behaviour. 
In our model we also assume that drivers do not respond to 
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animals on roads by taking evasive action, with this essentially 
being equivalent to assuming that the visibility of animals is zero. 
Incorporating visibility and the possibility of drivers taking evasive 
action would increase the probability of survival and would have 
an analogous effect to reducing traffic volumes. This effect would 
make it more likely for road density to have a larger impact that 
traffic volumes (Figure 3). Therefore, our results are also likely to 
be robust to variation in visibility across species. 

We quantified the relative impact on mortality from increases in 
traffic volume on existing roads and from increasing road density, 
but did not explicidy consider the spatial drivers of these patterns. 
The effect of each of these will depend upon the spatial locations of 
existing roads and the potential locations of any future roads 
relative to habitat. Where roads occur in close proximity to habitat 
or resources required by species, then this can result in the 
attraction of wildlife to roads, resulting in higher mortality rates 
[70,71]. Friar et al. [19] show that the road mortality hazard for 
elk (Cerous ekphus) is higher for roads associated with clearcuts than 
for roads independent of clearcuts. Therefore, the location of new 
roads is likely to be a critical factor in determining the actual 
impact on mortality. If new roads are necessary to accommodate 
increased traffic, then the impact on wildlife may be reduced by 
avoiding locations close to existing habitat. An alternative and 
potentially complementary strategy that has been proposed for 
reducing road impacts involves traffic calming in key habitat areas 
and redirecting traffic to surrounding roads [20,48]. For example, 
van Langevelde and Jaarsma [20] show that traffic calming results 
in improvements in roe deer (Capreolus capreolus) persistence by 
reducing fragmentation effects. The mechanism underlying the 
success of this strategy is similar to the mechanism that drives the 
relative success of accommodating increased vehicle numbers by 
raising capacity on existing roads versus building new roads; it 
reduces the number of high-volume roads that individual animals 
need to cross. However, the effect of traffic calming is also likely to 
depend critically on the spatial distribution of both habitat and the 
selection of areas for traffic calming. Understanding how the 
spatial location of new roads and calmed areas affect wildlife in 
developing principles for decision-making present important 
challenges for future research. 

We restricted our analysis to investigating some simple 
principles about how planners can reduce the impact of roads 
on wildlife by considering the relative effects of road density and 

References 

1. Roger E, Laffan SW, Ramp D (201 1) Road impacts a tipping point tor wildlife 
populations in threatened landscapes. Popul Ecol 53: 215-227. 

2. Trombulak SC, Frissell CA (2000) Review of ecological effects of roads on 
terrestrial and aquatic communities. Conscrv Biol 14: 18-30. 

3. Tigas LA, Van Vuren DH, Sauvajot RM (2002) Behavioral responses of bobcats 
and coyotes to habitat fragmentation and corridors in an urban environment. 
Biol Conscrv 108: 299-306. 

4. Goosem M (1997) Internal fragmentation: the effects of roads, highways, and 
powerline clearings on movements and mortality of rainforest vertebrates. In: 
Laurence WF, Bierregaard RO, editors. Tropical forest remnants Ecology, 
management, and conservation of fragmented communities. Chicago, USA: The 
University of Chicago Press, pp. 241—255. 

5. Dudaniec RY, Rhodes JR, Worthington Wilmer J, Lyons M, Lee KE, et al. 
(2013) Using multi-level models to identify drivers of landscape genetic structure 
among management areas. Mol Ecol 22: 3752-3765. 

6. Gibbs JP, Shrivcr WG (2002) Estimating the effects of road mortality on turtle 
populations. Conscrv Biol 16: 1647-1652. 

7. Wilkic D, Shaw E, Rotberg F, Morclli G, Auzcl P (2000) Roads, development, 
and conservation in the Congo basin. Conscrv Biol 14: 1614—1622. 

8. Lunney D (in press) Wildlife roadkill: illuminating and overcoming a blind spot 
in public perception. Pac Conscrv Biol. 

9. Hels T, Buehwald E (2001) The effect of road kills on amphibian populations. 
Biol Conscrv 99: 331-340. 

10. Ramp D, Ben-Ami D (2006) The effect of road-based fatalities on the viability of 
a peri-urban swamp wallaby population. J Wildl Manag 70: 1615-1624. 



traffic volume. However, these are not the only tools available for 
reducing road impacts. Numerous other strategies are regularly 
employed to help reduce the risk of wildlife mortality on existing 
roads, including speed reduction measures, fencing, road crossings, 
culverts under roads, overpasses, lighting, signs, road threshold 
treatments and wildlife reflectors, all with varying degrees of 
success [16]. In Port Stephens, a combination of advisory speed 
reductions zones, fencing and koala crossings have been proposed 
as part of a koala plan of management for the local government 
area [44] . However, the relative effectiveness of strategies that aim 
to reduce mortalities on existing roads compared to strategies 
focused on road network design remains uncertain. Addressing this 
would require not only an analysis of the effectiveness of each 
strategy, but also an estimation of the costs of implementation to 
identify the most cost-effective approach to reducing wildlife 
mortality [72,73]. An important area for future work is therefore 
to incorporate costs into models of the effectiveness of alternative 
strategies for mitigating road mortalities, so as to identify 
investment priorities that achieve mortality reductions for lowest 
cost. 

The impact of roads on wildlife populations arises from a range 
of complex spatial processes involving interactions between 
movement behaviours and the spatial pattern of habitat and 
roads. This study makes an important contribution to under- 
standing how best to accommodate future increases in vehicle 
numbers and presents a coherent approach for doing this. One of 
the key challenges now is developing ways to effectively integrate 
the results of studies such as this into strategic planning processes 
for infrastructure and wildlife management [74] . This is challeng- 
ing, but vital for biodiversity conservation. 

Acknowledgments 

Thank you to Steve Wilson (formerly of Port Stephens Council), to New 
South Wales Roads and Maritime Services for the road traffic data, and to 
Deidre de Villiers and the Queensland Department of Environment and 
Heritage Protection for the koala head to tail length data. 

Author Contributions 

Conceived and designed the experiments: JRR DLJC CAM. Performed 
the experiments: JRR. Analyzed the data: JRR. Contributed reagents/ 
materials/analysis tools: DL. Wrote the paper: JRR DLJC CAM. 



11. Fahrig L, Pedlar JH, Pope SE, Taylor PD, WegnerJF (1995) Effect of road 
traffic on amphibian density. Biol Conserv 73: 177—182. 

12. Jones ME (2000) Road upgrade, road mortality and remedial measures: impacts 
on a population of eastern quolls and Tasmanian devils. Wildlife Australia 27: 
289-296. 

13. Spellerbcrg IF (2002) Ecological effects of roads. Enfield, USA: Science 
Publishers. 

14. International Road Federation (201 1) World Road Statistics 201 1 - Data 2004 - 
2009. Geneva, Switzerland: International Road Federation. 

15. Fahrig L, Rytwinski T (2009) Effects of roads on animal abundance: an 
empirical review and synthesis. Ecol Soc 14. 

16. Glista DJ, DeVault TL, DeWoodyJA (2009) A review of mitigation measures for 
reducing wildlife mortality on roadways. Landsc Urban Plann 91: 1—7. 

17. Lesbarrercs D, Fahrig L (2012) Measures to reduce population fragmentation by 
roads: what has worked and how do we know? Trends Ecol Evol 27: 374-380. 

18. Taylor BD, Goldingay RL (2003) Cutting the carnage: wildlife usage of road 
culverts in north-eastern New South Wales. Wildl Res 30: 529-537. 

19. Frair JL, Merrill EH, Beyer HL, Morales JM (2008) Thresholds in landscape 
connectivity and mortality risks in response to growing road networks. J Appl 
Ecol 45: 1504-1513. 

20. van Langevelde F, Jaarsma CF (2009) Modeling the effect of traffic calming on 
local animal population persistence. Ecol Soc 14. 

2 1 . Jaeger JAG, Fahrig L, Ewald KC (2006) Does the configuration of road networks 
influence the degree to which roads affect wildlife populations? In: Irwin CL, 
Garrett P, McDcrmott NC, editors. Proceedings of the 2005 International 



PLOS ONE | www.plosone.org 



9 



March 2014 | Volume 9 | Issue 3 | e91093 



Designing Road Networks for Wildlife Conservation 



Conference on Ecology and Transportation. Raleigh, USA: Centre for 
Transportation and the Environment, North Carolina State University, pp. 
151-163. 

22. Dique DS, Thompson J, Prcece HJ, Penfold GC, de Villiers DL, et al. (2003) 
Koala mortality on roads in south-east Queensland: the koala speed-zone trial. 
Wildl Res 30: 419-426. 

23. Rhodes JR, Ng CF, de Villiers DL, Preece HJ, McAlpine CA, et al. (201 1) Using 
integrated population modelling to quantify the implications of multiple 
threatening processes for a rapidly declining population. Biol Conscrv 144: 
1081-1088. 

24. Natural Resource Management Ministerial Council (2009) National Koala 
Conservation and Management Strategy 2009-2014. Canberra, Australia: 
Commonwealth of Australia. 

25. Lunncy D, Close R, Crowther MS, Bryant J, Shannon I, et al. (2010) The koalas 
of Campbclltown, south-western Sydney: does their natural history foretell of an 
unnatural future? In: Lunney D, Hutchings P, Hochuli D, editors. The Natural 
History of Sydney. Mosman, Australia: Royal Zoological Society of NSW. 

26. DECC (2008) NSW (2008) Koala Recovery Plan. Sydney, Australia: 
Department of Environment and Climate Change. 

27. Rhodes JR, McAlpine CA, Lunney D, Possingham HP (2005) A spatially explicit 
habitat selection model incorporating home range behavior. Ecology 86: 1 199- 
1205. 

28. Hindcll MA, Handasyde KA, Lee AK (1985) Tree species selection by free- 
ranging koala populations in Victoria. Aust Wildl Res 12: 137-144. 

29. Callaghan J, McAlpine C, Thompson J, Mitchell D, Bowen M, et al. (2011) 
Ranking and mapping koala habitat quality for conservation planning on the 
basis of indirect evidence of tree-species use: a case study of Noosa Shire, south- 
eastern Queensland. Wildl Res 38: 89-102. 

30. White NA (1999) Ecology of the koala {Phascolarctos cinereus) in rural south-cast 
Queensland, Australia. Wildl Res 26: 731-744. 

31. Phillips S, Callaghan J, Thompson V (2000) The tree species preferences of 
koalas (Phascolarctos cinereus) inhabiting forest and woodland communities on 
Quaternary deposits in the Port Stephens area, New South Wales. Wildl Res 27: 
1-10. 

32. Hindell MA, Lee AK (1987) Habitat use and tree preferences of koalas in a 
mixed eucalypt forest. Aust Wildl Res 14: 349-360. 

33. Lunncy D, Phillips S, Callaghan J, Coburn D (1998) Determining the 
distribution of koala habitat across a shire as a basis for conservation: a case 
study from Port Stephens, New South Wales. Pae Conscrv Biol 4: 186—196. 

34. Lunncy D, Matthews A, Moon C, Ferricr S (2000) Incorporating habitat 
mapping into practical koala conservation on private lands. Conserv Biol 14: 
669-680. 

35. Cork SJ, Hume ID, Foley WJ (2000) Improving habitat models and their utility 
in koala conservation. Conserv Biol 14: 660-668. 

36. Mitchell P (1990) The home ranges and social activity of koalas - a quantitative 
analysis. In: Lee AK, Handasyde KA, Sanson CD, editors. Chipping Norton, 
Australia: Surrey Beatty & Sons. pp. 171-188. 

37. Ellis WAH, Mclzer A, Carrick FN, Hascgawa M (2002) Tree use, diet and home 
range of the koala [Phascolarctos cinereus) at Blair Athol, central Queensland. Wildl 
Res 29: 303-311. 

38. Seabrook L, McAlpine C, Baxter G, Rhodes J, Bradley A, et al. (20 1 1) Drought- 
driven change in wildlife distribution and numbers: a case study of koalas in 
south west Queensland. Wildl Res 38: 509-524. 

39. Adams-Hosking C, McAlpine C, Rhodes JR, Grantham HS, Moss PT (2012) 
Modelling changes in the distribution of the critical food resources of a specialist 
folivore in response to climate change. Divers Distrib 18: 847—860. 

40. Lunncy D, Grcsser S, O'Neil LE, Matthews A, Rhodes JR (2007) The impact of 
fire and dogs on koalas at Port Stephens, New South Wales, using population 
viability analysis. Pac Conserv Biol 13: 189—201. 

41. Rhodes JR, Wiegand T, McAlpine CA, Callaghan J, Lunney D, ct al. (2006) 
Modeling species' distributions to improve conservation in semiurban land- 
scapes: koala case study. Conserv Biol 20: 449^-59. 

42. Reed PC, Lunney D, Walker P (1990) A 1986-1987 survey of the koala 
Phascolarctos cinereus (Goldfuss) in New South Wales and an ecological 
interpretation of its distribution. In: Lee AK, Handasyde KA, Sanson GD, 
editors. Chipping Norton, Australia: Surrey Beatty & Sons. pp. 55-74. 

43. Knott T, Lunney D, Coburn D, Callaghan J (1998) An ecological history of 
koala habitat in Port Stephens Shire and the Lower Hunter on the central coast 
of New South Wales, 1801-1998. Pac Conserv Biol 4: 354-368. 

44. Port Stephens Council (2001) Port Stephens Council comprehensive koala plan 
of management. Part 1: the CKPoM. Raymond Terrace, Australia: Port 
Stephens Council. 



45. Melzer A (1995) Aspects of the ecology of the koala, Phascolarctos cinereus 
(Goldfuss, 1817), in the sub-humid woodlands of central Queensland. Brisbane, 
Australia: The University of Queensland. 

46. Pittman J (1993) Probability. New York: Springe r-Verlag. 

47. van Langeveldc F, Jaarsma CF (2004) Using traffic flow theory to model traffic 
mortality in mammals. Landsc Ecol 19: 895—907. 

48. van Langeveldc F, Jaarsma C (1997) Habitat fragmentation, the role of rural 
roads and their travcrsability. In: Canters K, editor. Netherlands: Ministry of 
Transport, Public Works and Water Management, pp. 171—182- 

49. Zar JH (1996) Biostatistical Analysis. Upper Saddle River, USA: Prentice-Hall. 

50. Wade PR (2000) Bayesian methods in conservation biology. Conscrv Biol 14: 
1308-1316. 

51. Stewart-Oaten A (1996) Goals in environmental monitoring. In: Schmitt RJ, 
Oscnburg CW, editors. Detecting ecological impacts: concepts and applications 
in coastal habitats. San Diego, USA: Academic Press, pp. 17-27. 

52. Worton BJ (1989) Kernel methods for estimating the utilization distribution in 
home-range studies. Ecology 70: 164-168. 

53. Hosmer DW, Lcmcshow S (2000) Applied logistic regression. New York, USA: 
John Wiley & Sons. 

54. Hosmer DW, Hosmer T, lcCessic S, Lcmcshow S (1997) A comparison of 
goodncss-of-fit tests for the logistic regression model. Stat Med 16: 965-980. 

55. Osius G, Rojek D (1992) Normal goodncss-of-fit tests for multinomial models 
with large degrees of freedom. Journal of the American Statistical Association 
87: 1145-1152. 

56. Slootcn E, Fletcher D, Taylor BL (2000) Accounting for uncertainty in risk 
assessment: case study of Hector's dolphin mortality due to gillnet entanglement. 
Conserv Biol 14: 1264-1270. 

57. Wade PR (2002) Bayesian population viability analysis. In: Bcissingcr SR, 
McCullough DR, editors. Population viability analysis. Chicago, USA: 
University of Chicago Press, pp. 213-238. 

58. Efron B, Tibshirani RJ (1993) An introduction to the bootstrap. New York, 
USA: Chapman & Hall. 

59. Caswell H (2001) Matrix population models. Sunderland, USA: Sinauer 
Associates. 

60. SushinskyJR, Rhodes J R, Possingham HP, Gill TK, Fuller RA (2013) How 
should we grow cities to minimize their biodiversity impacts? Global Change 
Biology 19: 401-410. 

61. McKinney ML (2002) Urbanization, biodiversity, and conservation. Bioscience 
52: 883-890. 

62. Forman RTT, Sperling D, Bissonctte JA, Clevenger AP, Cutshall CD, ct al. 
(2003) Road Ecology. Science and Solutions. Washington DC, USA: Island 
Press. 

63. Fields G, Hartgen D, Moore A, Poole RW, Jr. (2009) Relieving congestion by 
adding road capacity and tolling. International Journal of Sustainable 
Transportation 3: 360-372. 

64. Jaeger JAG, Bowman J, Brcnnan J, Fahrig L, Bert D, et al. (2005) Predicting 
when animal populations are at risk from roads: an interactive model of road 
avoidance behavior. Ecol Model 185: 329-348. 

65. Gunson KE, Mountrakis G, Quackcnbush LJ (2011) Spatial wildlife-vehicle 
collision models: a review of current work and its application to transportation 
mitigation projects. J Environ Manag 92: 1074-1082. 

66. Forman RTT, Alexander LE (1998) Roads and their major ecological effects. 
Annu Rev Ecol Syst 29: 207-231. 

67. Clarke GP, White PCL, Harris S (1998) Effects of roads on badger Meles meles 
populations in south-west England. Biol Conscrv 86: 1 17-124. 

68. Rondinini C, Doncaster CP (2002) Roads as barriers to movement for 
hedgehogs. Funct Ecol 16: 504-509. 

69. Bautista LM, Garcia JT, Calmaestra RG, Palacin C, Martin CA, et al. (2004) 
Effect of weekend road traffic on the use of space by raptors. Conserv Biol 18: 
726-732. 

70. Ramp D, Caldwell J, Edwards KA, Warton D, Croft DB (2005) Modelling of 
wildlife fatality hotspots along the snowy mountain highway in New South 
Wales, Australia. Biol Conserv 126: 474-490. 

71. Grosman PD, Jaeger JAG, Biron PM, Dussault C, OuelletJ-P (2011) Trade-off 
between road avoidance and attraction by roadside salt pools in moose: An 
agent-based model to assess measures for reducing moosc-vchiele collisions. Ecol 
Model 222: 1423-1435. 

72. Murdoch W, Polansky S, Wilson KA, Possingham HP, Kareiva P, ct al. (2007) 
Maximising return on investment in conservation. Biol Conscrv 139: 375-388. 

73. Polak T, Rhodes JR, Jones D, Possingham HP, (In press) Optimal planning for 
mitigating the impacts of roads on wildlife. J Appl Ecol. 

74. Dovers SR, Norton TW, Handmer JW (1996) Uncertainty, ecology, sustain- 
ability and policy. Biodivers Conserv 5: 1 143-1 167. 



PLOS ONE | www.plosone.org 



10 



March 2014 | Volume 9 | Issue 3 | e91093 



